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ABSTRACT 

The  wind  driven  surface  of  a  large  anechoic  tank  was  used 
to  study  forward  scattered  underwater  sound.   A  10  kHz  saw- 
tooth signal  was  used  to  drive  an  omnidirectional  source  to 
insonify  the  rough  surface.   Direct  and  reflected  path  sound 
signals  to  a  point  hydrophone  were  separated  by  appropriate 
gating  and  the  reflected  path  signal  was  analyzed  for  up  to 
2  minutes  at  10  and  20  msec,  intervals,  using  digital  FFT 
spectral  analysis.   Graphs  of  sound  pressure  level  vs.  time 
for  every  10  kHz  from  30  kHz  to  120  kHz  (Surface  roughnesses 
between  1.1  and  10.5)  are  shown.   Correlation  of  the  scattered 
sound  pressure  with  frequency  was  analyzed.   Evidence  is  pre- 
sented to  support  the  hypothesis  that,  for  large  roughnesses, 
the  instantaneous  scattered  sound  amplitudes  are  harmonic 
functions  of  sound  frequency;  that  is,  there  is  a  consistent 
frequency  separation  between  instantaneous  amplitude  maxima 
at  one  frequency  and  minima  at  another  frequency. 
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I.   INTRODUCTION 

For  a  number  of  years,  acousticians  have  studied  the 
fluctuations  of  underwater  acoustic  signals  when  scattered 
from  a  statistically  rough  interface,  typically  either  the 
surface  or  bottom.   This  effect  has  been  described  as  "build- 
ing and  fading"  of  the  signal;  an  analogy  has  frequently  been 
drawn  between  this  effect  and  the  twinkling  of  stars.   The 
mechanism  of  this  scattering  has  been  the  object  of  many 
theories;  References  1  through  5  refer  to  the  basic  develop- 
ments in  this  area.   The  theoretical  developments  have  common 
shortcomings :   in  order  to  solve  the  integral  equation  gov- 
erning the  scattered  sound  field,  it  is  necessary  to  make 
many  assumptions  involving  the  boundary  conditions,  the  value 
of  the  normal  derivative  on  the  surface,  the  spatial  statis- 
tics of  the  insonified  area;  most  analyses  deal  with  high  or 
low  frequency  limits;  shadowing  and  multiple  scattering  is 
usually  ignored  in  all  but  the  most  complex  and  sophisticated 
treatments.   These  comments  concerning  the  shortcomings  of 
strictly  theoretical  approaches  to  the  scattering  problem  have 
been  made  in  order  that  the  need  for  supplemental  experimental 
research  be  noted.   A  mixture  of  theoretical  insight  and  con- 
trolled experimental  procedure  is  sometimes  more  fruitful 
than  purely  theoretical  studies  with  arbitrary  assumptions. 
Such  is  the  philosophy  adopted  for  this  research. 
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A.  OBJECTIVES 

The  objective  of  this  research  is  to  investigate  the 
effects  that  variation  of  frequency  has  on  the  sound  scattered 
from  a  statistically  rough,  fluctuating,  wind  driven  water 
surface,  and  to  determine  if  a  predictable  relationship  exists 
between  the  intensities  of  various  harmonics  of  a  given  sound 
source  when  scattered  from  such  a  surface.   The  basic  con- 
cepts of  this  experimental  effort  are: 

1.  Insonification  of  the  fluctuating  surface  with  many 
frequencies  simultaneously. 

2.  Isolation  of  surface  reflected  sound  from  other  paths 
through  the  use  of  gated  signals. 

B.  INTERMEDIATE  GOALS 

The  following  intermediate  goals  are  presented  as  logical 
steps  in  attainment  of  the  objective  described  above: 

1.  Perform  preliminary  studies  of  the  wind  driven  surface, 
and  monof requency  sound  scattered  from  such  a  surface. 

2.  Develop  the  electronic  circuitry  to  transmit  a  wide 
range  of  frequencies  simultaneously,  receive  and  process 
these  signals. 

3.  Develop  and  test  an  acoustic  source  capable  of  trans- 
mitting over  a  wide  range  of  frequencies,  with  the  stip- 
ulation that  frequencies  of  interest  be  far  enough  above 
the  'noise'  to  permit  meaningful  signal  processing. 

4.  Analyze  the  processed  data  to  see  if  the  stated  objec- 
tive can  be  met. 
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II.   RESEARCH  FACILITIES 

A.  OCEAN  ACOUSTIC  WAVE  FACILITY  (OAWF) 

The  Ocean  Acoustic  Wave  Facility  is  a  combination  of  a 
wind  generated  wave  tunnel  and  an  anechoic  acoustic  tank. 
The  wind  tunnel  is  approximately  17  meters  long,  1.2  meters 
wide  and  1.2  meters  deep.   By  varying  the  number  of  fans  and/ 
or  the  distance  between  the  water  surface  and  the  top  of  the 
tunnel,  a  wide  range  of  ocean  surfaces  can  be  obtained. 
Figures  1  and  2  (from  Tourville)  show  the  essential  features 
of  the  OAWF. 

The  wave  generation  tunnel  empties  into  a  water  filled 
tank  with  dimensions  of  3  meters  on  each  side.   The  wind  gen- 
erated waves  are  dissipated  on  a  constructed  beach  to  prevent 
reflected  waves  from  interfering  with  wind  generated  waves. 
The  tank  is  constructed  of  4x4  redwood  posts,  placed  with 
corners  facing  outward.   The  use  of  redwood  provides  sound 
absorption,  as  shown  in  Reference  8;  however,  scattering  of 
sound  from  these  surfaces  did  present  some  experimental  diffi- 
culties, as  pointed  out  in  a  later  section.   Equipment  was 
positioned  within  the  tanks  by  use  of  movable  racks  as  shown 
in  Figure  3.   The  OAWF  models  the  real  ocean  on  a  scale  on  the 
order  of  150  to  1  up  to  sea  state  3. 

B.  DATA  ACQUISITION  AND  PROCESSING  CAPABILITIES 

Data  acquisition  and  processing  were  accomplished  by 
utilizing  a  digital  computer  system  composed  of  three  primary 
components  each  of  which  is  interfaced  to  provide  high  speed 
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analog  to  digital  conversion,  digital  processing,  and  data 
printout  and  recording.   The  design  was  developed  by  the 
Special  Projects  Section  of  the  Naval  Air  Development  Center 
in  conjunction  with  Pinkerton  Computer  Consultants,  Inc.,  of 
Warminster,  Pennsylvania.   The  three  components  are  pictured 
in  Figure  4  and  described  below. 

1.  Interdata  Model  70  Computer 

This  minicomputer  is  a  digital  design  that  is  FORTRAN 
programmable  with  a  64  thousand  byte  memory.   In  addition  to 
actual  core  memory,  data  that  have  been  stored  on  digital 
cassettes  can  be  read  into  the  computer  for  processing. 

2.  Phoenix  Analog  to  Digital  Converters,  Model  ADC  712 
Two  ADC  712  analog  to  digital  converters  may  be  used 

separately  or  simultaneously.   Each  converter  is  a  high  speed 
device  capable  of  encoding  +  10  volt  input  signals  in  digital 
form  and  providing  an  accuracy  of  0.005  volts  in  20  volts. 
The  maximum  sampling  rate  is  320,000  samples  per  second;  this 
rate  was  routinely  used. 

3.  Texas  Instruments  Silent  Electronic  Data  Terminal 
Model  73  3 

The  TI  733  consists  of  a  keyboard  used  as  a  program- 
ming input/output  control  device,  a  printer,  and  a  playback/ 
record  section  used  in  conjunction  with  digital  cassettes. 
The  overall  system  facilitates  rapid,  accurate  processing  of 
any  desired  type  of  analog  electrical  signal  and  was  used 
mainly  for  statistical  and  frequency  domain  analysis  using 
standard  FFT  algorithms.   Further  descriptions  are  contained 
in  References  6  and  7. 
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Figure  3.  Computer  System 


Figure  4.  View  of  OAWF 
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4.   IBM-360  Interface 

In  addition  to  the  digital  system  just  described, 
another  processing  capability  consisted  of  the  TI  733  termi- 
nal and  taped  digital  data  interfaced  via  phone  link  to  the 
IBM-360/67  CP/CMS  system.   This  capability  was  used  to  pro- 
vide data  rearrangement,  card  punch  capability  and  graphical 
output  (with  the  aid  of  a  CAL-COMP  plotter) . 


C.   STANDARD  EQUIPMENT  LIST 

Much  of  the  equipment  used  was  standard  off  the  shelf 
scientific  measurement  equipment;  a  list  of  standard  equip- 
ment referred  to  throughout  the  text  is  given  below: 


Abbreviation 
Scope 

HP  465A  Amp 

HP  4  66A  Amp 

HP  467A  Amp 

HP  218 

HP  219 

PAR  113  Amp 

KH  3322  Filter 
KH  3342  Filter 
KH  3350  Filter 
Freq.  Counter 


Full  Description 

Tectronix  Model  54 5B  4-channel 
oscilloscope 

Hewlett-Packard  Voltage  Amplifier 
Model  465A 

Hewlett-Packard  Voltage  Amplifier 
Model  466A 

Hewlett-Packard  Power  Amplifier 
Model  46  7A 

Hewlett-Packard  Model  218  Digital 
Delay  Generator 

Hewlett-Packard  Model  219  Dual 
Pulse  Unit 

Princeton  Applied  Research  Amplifier 
Model  113 

Krohn-Hite  Model  3322  Freq.  Filter 

Krohn-Hite  Model  3  34  2  Freq.  Filter 

Krohn-Hite  Model  3  350  Freq.  Filter 

Hewlett-Packard  5223L  Electronic 
Frequency  Counter 
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DMM 
GR-119  2B 

GR-1217C 

GR-139  6 

Impedance  Bridge 

Wavetek  114 
Wavetek  136 
Wavetek  144 


Fluke  Model  8000A  Digital  Multi  Meter 

General  Radio  Electronic  Frequency 
and  Period  Counter 

General  Radio  Model  1217C  Unit  Pulse 
Generator 

General  Radio  Model  139 6A  Tone 
Burst  Generator 

General  Radio  Model  16  50A  Impedance 
Bridge 

Wavetek  Model  114  Function  Generator 

Wavetek  Model  136  Function  Generator 

Wavetek  Model  144  Function  Generator 


Sound  sources,  receivers,  wave  probes  and  other  miscellaneous 
equipment  will  be  described  in  the  text,  as  necessary. 
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III.   EXPERIMENTAL  PROCEDURES 

The  experimental  procedures  followed  in  the  course  of 
this  research  can  be  grouped  into  four  general  categories: 

(1)  Calibration  of  the  Surface  Wave  Probe  and  determination 
of  wave  spectra,  wave  heights  and  surface  roughnesses; 

(2)  Continuous  wave  sound  scatter  measurements;  (3)  Sound 
Source  testing  and  selection;  (4)  Surface  Reflected  Pulsed 
Sound  Wave  analysis. 

A.   WAVE  PROBE  CALIBRATION,  WAVE  SPECTRA  AND  WAVE  HEIGHT 

MEASUREMENTS 

In  two  previous  works  [Refs.  6  and  7],  wave  height  fre- 
quency spectra  of  wind  generated  water  waves  were  determined. 
Confirmation  of  these  previous  studies  was  considered  to  be 
an  excellent  starting  point  of  this  current  research  and 
would  afford  valuable  insight  into  the  use  of  the  specialized 
measuring  equipment  and  into  the  use  and  programming  of  the 
digital  minicomputer.   It  was  intended  to  extend  these 
investigations  over  those  previously  performed  by  studying 
waves  generated  by  all  combinations  of  fans  and  also  to 
develop  a  program  to  give  actual  RMS  values  of  wave  height. 

1.   Wave  Probe  Calibration 

The  wave  height  probe  is  a  linear  capacitative  probe 
developed  by  Stanford  University  and  consists  basically  of  a 
thin  wire  whose  total  capacitance  is  determined  by  the  water 
depth  to  which  the  probe  is  immersed.   The  probe  is  connected 
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to  the  impedance  bridge  and  the  remainder  of  the  circuitry 
as  shown  in  Figure  5.   The  impedance  bridge  used  had  a  con- 
venient feature  which  provided  a  1  kHz  output  signal  modulated 
by  the  value  of  the  input  capacitance.   As  the  probe  depth  is 
varied  the  final  signal  (which  is  pure  D.C.  for  a  smooth 
surface)  provides  a  linear  output  proportional  to  the  probe 
depth.   The  probe  was  mounted  on  a  rigid  frame,  but  was  able 
to  be  varied  in  depth  by  means  of  a  rotating  knob-screw 
arrangement.   A  graduated  depth  scale  was  attached  to  the  frame 
with  an  accuracy  of  about  0.02  cm.   The  probe  was  then  cali- 
brated as  shown  in  Figure  6.   Initial  attempts  at  calibration 
were  not  successful  due  to  lack  of  repeatability.   This  was 
corrected  by  slight  modifications  to  the  filtering  network 
and  the  final  calibration  curve  was  satisfactory. 
2 .   Wave  Spectra  Measurements 

Upon  completion  of  the  wave  height  probe  calibration, 
measurement  of  wave  height  spectra  was  undertaken  using  the 
circuitry  of  Figure  5  as  shown  previously.   The  time  varying 
signal  proportional  to  instantaneous  wave  height  was  input 
into  one  A/D  channel  of  the  digital  computer  for  analysis 
using  the  SKIP-1C  program  (all  programs  used  are  described 
briefly  in  Appendix  A) .   The  computer  A/D  sampling  rate  was 

g 

set  at  20  Hz;  50  runs  of  256  (2  )  samples  each  were  taken  at 
this  sample  rate.   Previous  spectral  measurements  had  indi- 
cated that  most  of  the  spectral  components  were  below  4  Hz 
so  that  the  Nyquist  criterion  was  satisfied  in  all  cases. 
Spectral  resolution  for  this  particular  data  gathering 
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Figure  5.   Wave  Height  Circuitry 
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arrangement  is  0.07  8  Hz.   For  each  of  the  50  data  sets, 
Fast  Fourier  Transform  (FFT)  analysis  was  performed  and  the 
output  spectral  density  for  each  run  was  computed,  normalized 
and  averaged.   The  output  spectral  densities  of  surfaces  for 
all  fan  combinations  of  one  to  five  fans  were  computed  and 
plotted.   The  graphs  of  output  power  spectral  densities  vs. 
frequencies  are  shown  in  Figure  7.   This  set  of  curves  shows 
the  expected  spectral  shifts  as  the  wave  heights  and  wave 
periods  are  increased  by  the  use  of  more  fans.   The  general 
results  were  in  agreement  with  Perkins  and  Tourville,  hence 
further  data  processing  simply  for  the  purpose  of  curve  smooth- 
ing was  not  performed.   As  a  further  check  on  the  computer 
accuracy  and  program,  signals  were  also  input  to  a  PAR  Fourier 
Analyzer  and  the  spectral  output  vs.  frequency  were  plotted 
on  an  XY  recorder.   Spectra  produced  by  this  alternate  means 
were  consistent  with  those  analyzed  by  digital  means.   A  typ- 
ical example  of  this  alternate  spectral  curve  is  shown  in 
Figure  8. 

3.   Average  Wave  Height  Determination 

After  the  wave  height  probe  calibration  curve  was 
obtained,  the  next  task  to  be  undertaken  was  accurate  deter- 
mination of  average  wave  heights  for  all  fan  combinations. 
This  was  accomplished  by  using  the  circuitry  of  Figure  5 
once  again,  in  conjunction  with  a  computer  program  entitled 
CAL  01.   This  program  sampled  the  wave  height  output  voltage, 
performed  root-mean-square  averaging  of  the  total  number  of 
A/D  data  points,  and  compared  this  with  a  reference  sin  wave 
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voltage  set  exactly  to  1.000  volts.   (This  comparison  elimi- 
nated the  need  for  A/D  voltage  calibration  each  time  the 
program  was  used.)   This  program  was  extensively  tested  with 
signals  of  known  RMS  values  to  ensure  accuracy  prior  to 
experimentally  determining  the  wave  heights.   Once  the  pro- 
gram was  known  to  be  working  properly,  RMS  voltage  values 
for  all  fan  combinations  were  determined,  and  converted  to 
wave  heights  using  the  linear  calibration  curve.   These  values 
are  shown  in  Table  I.   Once  the  average  wave  heights  were 

Table  I 
RMS  Wave  Heights  vs.  Number  of  Fans 

No.  of  Fans           AC  (RMS)  Volts  RMS  Wave  Height  (Cm) 

1  0.01475  0.447 

2  0.0180  0.546 

3  0.0244  0.740 

4  0.0343  1.041 

5  0.0416  1.262 

obtained  for  each  fan  combination,  it  was  then  a  simple  matter 
to  compute  surface  roughness  values  using  the  definition: 

h                2ttq    ._ 
g   =  — r (2  cos  9) 

where  g   is  defined  as  the  roughness  parameter,  a  is  the  RMS 
wave  height,  X   is  the  wavelength  of  the  sound  source  and  9 
is  the  angle  of  incidence.   The  roughness  values  were  com- 
puted using  the  geometry  of  Figure  21;  values  are  contained 
in  Table  II. 
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Table  II 
Roughness  Parameters 


Freq  (kHz) 

A  Cm) 

Roughness  -yg" 

10 

0.15 

1  Fan 

2  Fans 

3  Fans 

4  Fans 

0.37 

0.45 

0.62 

0.86 

20 

0.075 

0.74 

0.91 

1.23 

1.72 

30 

0.050 

1.11 

1.36 

1.85 

2.60 

40 

0.0375 

1.49 

1.81 

2.46 

3.51 

50 

0.030 

1.86 

2.27 

3.08 

4.33 

60 

0.0250 

2.23 

2.72 

3.69 

5.20 

70 

0.0214 

2.60 

3.18 

4.31 

6.06 

80 

0.01875 

2.97 

3.63 

4.92 

6.92 

90 

0.0166 

3.35 

4.08 

5.54 

7.79 

100 

0.0150 

3.72 

4.53 

6.15 

8.7 

110 

0.0136 

4.09 

4.99 

6.77 

9.52 

120 

0.0125 

4.46 

5.45 

7.38 

10.38 

B.   CONTINUOUS  WAVE  SOUND  SCATTER  MEASUREMENTS 

It  was  decided  to  initially  analyze  the  scattered  sound 
spectra  of  a  continuous  20  kHz  source.   The  source  and  re- 
ceiver setup  and  associated  circuitry  are  shown  in  Figure  9 
and  are  similar  to  that  of  Perkins  so  that  comparison  of  re- 
sults could  be  made.   The  source  used  was  a  2"  ceramic 
sphere;  the  receiver  was  a  standard  LC-10  hydrophone.   The 
receiver  was  positioned  in  a  Lloyds1  mirror  maximum  signal 
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Figure  9. Sound  Processing  Circuitry 
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location.   The  demodulated  scattered  sound  envelope  was  input 
into  the  digital  computer  and  analyzed  with  the  SKIP-1C  pro- 
gram.  The  normalized  scattered  sound  spectral  densities 
for  1  to  5  fans  are  shown  in  Figures  10-14.   Each  spectral 
density  curve  is  an  average  of  50  FFT's,  each  consisting  of 
256  sample  points,  sampled  at  a  rate  of  128  samples/sec. 
Once  again  results  were  consistent  with  previous  experimen- 
ters' results.   Oscilloscope  photos  of  typical  modulated  and 
demodulated  sound  signals  are  shown  in  Figures  15  and  16. 

C.   SOUND  SOURCE  SELECTION  AND  TESTING 

As  discussed  in  the  introductory  chapter,  in  order  to 
insonify  an  area  simultaneously  with  a  wide  range  of  frequen- 
cies, it  was  necessary  to  test  various  sources  to  see  which 
could  provide  the  maximum  transmitting  frequency  response  in 
the  range  of  interest,  nominally  5  to  160  kHz. 

1.   Transducer  Construction  and  Testing 

A  significant  number  of  potential  transducers  were 
constructed  and  tested;  for  the  sake  of  brevity  only  a  brief 
description  of  each  will  be  given  along  with  the  reason  for 
rejection. 

a.   Two-inch  ceramic  sphere,  hollow  interior:   this 
transducer  was  utilized  in  the  initial  sound  scattering 
studies  at  20  kHz.   The  transducer  did  not  have  a  significant 
output  below  20  kHz  and  significant  'ringing'  occurred  at 
50  kHz  which  would  have  made  the  sphere  unsuitable  for  short 
separation  gated  pulses. 
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FIGURE  13.  Surface  Modulated  Sound  Spectrum, 4  Fans. 
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Figure  15.  Surface  Modulated  Sound  Signal 


Figure  16.  Surface  Demodulated  Sound  Signal^ 
(Figure  16  is  not  the  demodulation 
of  Figure  15, but  was  taken  at  a 
different  time  and  is  a  typical 
demodulated  signal.) 
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b.  One-inch  ceramic  sphere,  hollow  interior:   same 
comments  as  in  a.,  extensive  ringing  around  75  kHz. 

c.  One-inch  ceramic  sphere,  oil  filled:   same  prob- 
lems as  previously  discussed.   Oil  was  introduced  into  the 
cavity  to  damp  out  the  ringing  phenomenon  but  it  only  shifted 
the  resonance  values. 

d.  One-inch  ceramic  sphere,  oil  filled  with  steel 
wool  damping.   No  major  improvement  noted  here. 

e.  One-inch  ceramic  sphere,  packed  with  heavy 
powdered  metal.   Same  comments  as  c.  apply. 

f.  Four  cylindrical  mylar  transducers  were  con- 
structed, 1"  and  2"  in  diameter,  10"  and  20"  in  length  respec- 
tively.  The  transducers  had  an  aluminum  core,  thin  mylar 
covering,  with  the  entire  assembly  coated  with  neoprene  water- 
proofing compound.   When  tested,  all  four  of  these  transducers 
exhibited  too  low  a  signal  level  to  be  useful. 

g.  LC-10  ceramic  hydrophone.   This  hydrophone  ex- 
hibited satisfactory  signal  level  over  the  entire  frequency 
range;  due  to  its  small  size  and  relative  omnidirectional 
features,  this  standard  receiving  hydrophone  was  utilized  as 
an  acceptable  source  throughout  the  remainder  of  the  experi- 
ment.  Actual  characteristics  of  this  transducer  will  be 
shown  and  discussed  in  subsequent  sections. 

D.   SURFACE  REFLECTED  PULSED  WAVE  ANALYSIS 

This  section  will  detail  the  experimental  procedure  used 
to  analyze  the  surface  reflected  sound  as  a  function  of 
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frequency  and  time.  Prior  to  studying  the  behavior  of  sound 
in  the  anechoic  tank,  a  study  of  the  sawtooth  oscillator  was 
initiated  as  described  below. 

1.   Wavetek  14  4  Sawtooth  Waveform  Analysis 

As  previously  mentioned,  a  Wavetek  144  function  gen- 
erator was  modified  to  give  an  asymmetrical  sawtooth  wave- 
with  a  slope  ratio  of  44:1.   This  asymmetrical  output  was 
coupled  directly  into  the  A/D  input  channel  of  the  digital 
computer.   The  THCDB  program  was  utilized  for  Fourier  analyses 
of  this  input  signal.   Each  entire  analysis  of  the  10  kHz 
sawtooth  waveform  consisted  of  4  9  sample  blocks,  each  block 
consisting  of  256  sample  points  taken  at  a  sample  rate  of 
320  kHz,  giving  an  available  frequency  resolution  of  1250  Hz. 
The  sawtooth  waveform  was  analyzed  in  both  a  continuous  and 
gated  mode.   Examination  of  the  output  data  (mean  spectral 
density  values,  and  standard  deviation  vs.  5  kHz  frequency 
increments)  revealed  an  unacceptably  high  value  in  the  higher 
frequency  ranges .   Values  of  standard  deviation  in  dB  re 
1  ybar  for  frequencies  over  100  kHz  were  greater  than  3  dB 
in  all  cases.   In  order  to  make  subsequent  experiments  mean- 
ingful, a  method  had  to  be  found  to  decrease  the  standard 
deviation  to  less  than  1  dB  for  all  frequencies  to  be 
examined. 

The  first  study  performed  was  to  compare  the  Wavetek 
144  spectral  analyses  with  similar  results  from  other  signal 
generators.   As  can  be  seen  in  Table  III,  three  different 
models  of  function  generators  gave  similar  results,  showing 
the  increased  instability  of  harmonics  with  frequency. 
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Various  methods  were  tried  in  order  to  overcome  this 
difficulty;  amplitude  variation  was  investigated  but  had  no 
significant  effect.   The  ultimate  solution  to  this  problem 
proved  to  be  band  pass  filtering  of  the  output  of  the  Wavetek 
prior  to  analysis.   An  extensive  study  of  the  effects  of  var- 
iation of  filter  limits  was  performed;  results  showed  that 
filtering  was  effective  in  reducing  the  standard  deviation 
of  the  output  signal  to  within  0.5  dB  or  less.   Further 
studies  of  this  nature  were  performed  with  the  direct  and 
reflected  sound  paths  in  the  wave  tank  and  will  be  subse- 
quently discussed. 

2.   Circuit  Arrangement  for  Sound  Analysis 

The  basic  circuit  arrangement  for  all  future  dis- 
cussions on  sound  analysis  is  shown  in  Figure  17.   A  brief 
discussion  of  the  basic  function  of  each  component  is  neces- 
sary to  understand  the  processing  of  the  various  sound  sig- 
nals.  The  Wavetek  144  function  generator  controls  the  basic 
repetition  frequency  of  the  entire  transmissions  circuit. 
The  GR  1217  Pulse  Generator  controls  the  actual  length  of  the 
sawtooth  waveform,  and  the  positive  going  portion  of  its 
pulse  initiates  the  timing  sequence  of  the  HP  218/219.   The 
HP  218/219  Digital  Delay  Generator,  when  operated  in  the 
'common'  mode  has  as  its  output  two  pulses  of  variable  ampli- 
tude and  width,  each  of  which  can  be  set  to  a  pre-adjusted 
time  delay  in  microseconds.   This  component  also  has  a  syn- 
chronizing output  for  use  in  oscilloscope  'synch'  input.   The 
time  delayed  pulses  are  used  as  the  timing  or  gating  signals 
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for  the  GR-1394  Tone  Burst  generator.   This  component,  in 
turn,  allows  the  passage  of  the  320  kHz  sampling  frequency 
to  the  computer  only  as  the  gating  pulses  permit. 

The  GR-1217  also  controls  the  sawtooth  waveform 
going  to  the  LC-10  source.   The  receiving  circuit  consists 
of  an  LC-10  hydrophone,  80  dB  of  amplification,  and  a  band- 
pass filter.   The  output  signal  of  the  filter  goes  to  the 
A/D  input  channel  for  computer  analysis.   The  series  of 
oscilloscope  photographs  shown  in  Figures  18-20  present  the 
various  waveforms  and  their  interrelation.   In  examination 
of  these  figures,  the  appearance  of  the  waveforms  should  be 
noted.   In  Figures  19  and  20,  the  direct  path  and  the  surface 
reflected  path  can  easily  be  identified;  however  in  both 
cases,  additional  scattering  is  present  following  the  surface 
reflected  wave,  consisting  of  a  mixture  of  acoustic  energy 
scattered  from  the  side  and  side-surface  of  the  wave  tank. 
It  was  the  presence  of  this  additional  scattered  wave  which 
determined  the  minimum  separation  (about  7  msec.)  between 
successive  sound  pulses  and  also  determined  the  maximum  length 
the  surface  reflected  pulse  could  attain  (about  0.6  msec.) 
without  overlap  with  the  scattering. 
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CONTROL  PULSE, GR- 12 17  OUTPUT 


WAVETEK  144  GATED  10  KHZ  OUT 


HP- 218/219  DELAYED  PULSES 


GATED  3  20KHZ  COMPUTER  SAMPLING 
SIGNAL.  GR-1396  OUTPUT. 


Figure  18 


WAVETEK  144  GATED  OUTPUT 


Figure  19 


DIRECT  AND  REFLECTED  SAWTOOTH 
SOUND  SIGNALS, FOLLOWED  BY 
SCATTERING. 

COMPUTER  SAMPLING  SIGNAL 
LOCATED  DIRECTLY  BELOW 
SURFACE  REFLECTED  SOUND 
SIGNAL. 


SERIES  OF  DIRECT, REFLECTED, 
AND  SCATTERED  SOUND  PULSES 
SHOWN  WITH  A  20  MSEC.  SEPARA- 
TION. 


Figure  20 
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IV.   EXPERIMENTAL  RESULTS  AND  ANALYSIS 

A.   STATIC  ACOUSTIC  WAVEFORM  INVESTIGATION 

Upon  completion  of  the  testing  and  spectral  analysis  of 
the  direct  output  of  the  Wavetek  144  oscillator,  testing  of 
the  direct  path  and  smooth  surface  reflected  acoustic  sig- 
nals was  initiated.   Although  some  variation  of  geometry  was 
initially  utilized,  the  geometry  of  Figure  21  was  used  for 
all  analysis  unless  otherwise  noted.   The  circuit  arrange- 
ment of  Figure  17,  shown  in  the  previous  section,  applies. 

A  study  of  the  spectral  components  of  a  10  kHz  sawtooth 
waveform  transmitted  via  the  two  different  paths  of  interest 
was  conducted.   Although  filtering  of  the  signal  was  per- 
formed, the  old  problem  of  high  standard  deviation  among  the 
components  was  noted  but  could  not  be  reduced  to  acceptable 
levels;  many  of  the  standard  deviations  were  on  the  order  of 
2  dB  and  this  was  felt  to  be  unacceptable  for  the  supposedly 
stable  direct  and  reflected  path  signals.   Further  study  re- 
vealed that  both  the  PAR- 113  and  KH-3550  had  voltage  limita- 
tions on  their  output,  which  if  exceeded,  resulted  in  clip- 
ping of  the  output.   The  PAR- 113  started  to  clip  at  ±  6  volts 
and  the  KH-3550  at  an  output  of  i  8  volts.   By  reducing  the 
magnitude  of  the  input  signal,  the  amplified  received  sound 
could  be  kept  within  these  bounds;  when  the  signals  were  kept 
within  these  limits  standard  deviations  were  greatly  reduced 
and  were  acceptable  (Table  IV) .   Due  to  the  inherent  limita- 
tions of  the  PAR- 113,  it  was  replaced  by  a  HP-465  amplifier 
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Figure  21.  Source  and  Receiver  Geometry 
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which  had  no  limitations  near  the  range  of  voltage  values 
encountered.   Particular  care  had  to  be  taken  when  analyzing 
the  fluctuating  signals  of  the  rough  surface  to  ensure  the 
clipping  point  of  the  KH-3550  was  not  approached.   During 
this  phase  of  testing,  the  filter  limits  of  20-130  kHz  were 
•found  to  produce  the  most  stable  signals  with  the  least 
overall  standard  deviation  of  the  direct  and  reflected 
waveforms. 

The  spectral  output  of  the  direct  and  surface  reflec- 
ted sound  paths,  as  analyzed  by  the  THCDB  program  are  shown 
in  Table  IV  below.   These  values  are  representative  for  the 
smooth  surface  case.   Typical  noise  values,  obtained  by  dis- 
connecting the  source,  and  analyzing  the  received  signal  with 
the  THCDB  program  are  also  shown.   The  results  of  this  table 
are  shown  graphically  in  Figure  22.   Each  mean  value  and 
standard  deviation  were  computed  over  a  50  block  average,  with 
each  block  consisting  of  128  sample  points  sampled  at  320  kHz. 
The  slightly  higher  standard  deviation  for  the  reflected  path 
at  120  kHz  is  thought  to  be  due  to  a  minor  surface  rippling 
due  to  slight  air  movement  over  the  wave  tank.   Note  that  most 
signal  to  noise  ratios  for  the  reflected  wave  are  on  the  order 
of  30  dB  or  higher. 

B.   ROUGH  SURFACE  SOUND  SCATTER  MEASUREMENT 

Surface  sound  scatter  measurements  were  now  attempted  for 
all  combinations  of  1  to  5  fans.   As  before,  the  circuit 
arrangement  of  Figure  17  was  used.   Gated  sound  pulses  were 
set  to  a  20  msec,  interval  between  pulses;  the  A/D  sampling 
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pulses  were  positioned  in  time  to  be  directly  coincident  with 
reception  of  the  surface  scattered  sound  pulses.   Each  of  the 
50  blocks  consisted  of  128  A/D  samples  taken  at  a  rate  of 
320  kHz;  thus  each  sample  block  took  0.4  msec.   During  this 
small  time  interval,  it  was  assumed  that  the  rough  wave  sur- 
face is  frozen  in  space  over  the  sampling  period.   This  assump- 
tion is  justified  based  on  examination  of  the  wave  frequency- 
spectra  previously  shown  where  the  peak  energy  is  at  a  period 
of  approximately  300  msec. 

The  output  sound  pressure  levels  vs.  frequency  for  each  of 
the  time  separated  reflected  pulses  was  computed  using  the 
THCDB  program,  and  stored  on  a  separate  digital  tape  cassette 
for  each  fan  combination.   This  tape  was,  in  turn,  processed 
by  another  program,  SORTF,  which  sorted  the  SPL  vs.  time  by 
frequency.   At  this  point,  for  each  5  kHz  increment  of  frequency 
there  existed  a  set  of  SPL  data  over  50  points,  each  point 
separated  by  20  msec,  in  time,  encompassing  a  total  time  span 
of  1  sec.   The  data  of  the  SORTF  program  was  also  stored  on 
digital  tape  cassettes.   Utilizing  a  phone  link  with  the  IBM 
360  computer  complex,  the  SORTF  data  was  card  punched  and  sub- 
sequently plotted  using  the  IBM  360  plotting  routines  and  the 
CAL  COMP  plotter. 

1.   Initial  Results 

Typical  results  from  the  previously  described  process 
are  shown  in  Figures  23  and  24.   The  particular  curves  shown 
are  for  relatively  low  surface  roughness  values.   Study  of  the 
whole  series  of  curves  obtained  for  the  1-5  fan  cases  indicated 
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the  existence  of  a  possible  problem.   The  pseudo  periodicity 
of  the  SPL  curves  was  on  the  order  of  12-15  Hz;  this  was 
approximately   three  times  that  expected  for  the  lower  rough- 
ness values .   A  complete  review  of  experimental  and  computa- 
tional procedures  was  initiated  to  seek  the  answer  to  this 
problem.   A  reexamination  of  the  THCDB  program  showed  that 
after  the  A/D  data  for  one  sample  had  been  stored,  Fourier 
analysis  and  all  remaining  computational  processes  were  being 
performed  prior  to  the  storage  of  the  next  A/D  data  set.   In 
effect,  what  resulted  was  a  skipping  of  the  timed  computer 
sampling  orders  until  the  computational  steps  were  performed, 
which  resulted  in  timed  samples  being  taken  only  every  60  to 
80  msec,  vice  the  20  msec,  separation  which  was  expected.   As 
a  result  of  this  finding,  it  was  decided  to  modify  the  exist- 
ing  programs  to  eliminate  the  problem. 
2 .   Program  Modifications 

In  order  to  get  data  at  the  precise  time  desired,  a 
number  of  program  modifications  were  made  as  described  below: 

a.   Raw  A/D  data  were  taken  and  stored  on  a  digital 
cassette  using  an  existing  program,  PERKlA.   In  order  to  ensure 
the  data  were  being  taken  at  the  exact  time  desired,  testing 
was  performed  on  10  kHz  triangle  wave  with  A/D  data  being 
taken  with  a  5  msec,  separation.   The  resulting  A/D  data  were 
manually  plotted  and  fitted  exactly  to  the  predicted  triangle 
form.   The  PERKlA  program  was  subsequently  modified  to  be  able 
to  take  up  to  100  data  blocks  of  12  8  points  each.   The  modi- 
fied program  was  called  PRKMOD. 
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b.  The  THCDB  program  was  modified  to  be  able  to 
accept  input  data  from  a  cassette  and  to  be  able  to  compute 
and  store  on  tape  the  output  spectral  analyses  of  up  to  100 
blocks  in  time.   The  modified  program  was  tested  on  a  static 
reflected  wave  from  the  smooth  surface.   This  was  compared 
with  the  original  program  and  nearly  identical  results  were 
obtained;  the  modified  program  was  labeled  THCDB (1). 

c.  The  SORTF  program  was  modified  to  accept  up  to 
100  data  blocks. 

3.   Final  Experimental  Data 

With  the  program  modifications  described  above  final 
experimental  data  taking  was  commenced.   The  first  set  of  data 
taken  was  for  fan  combinations  of  1  to  5,  pulse  and  sample 
time  separation  of  2  0  msec,  for  a  total  duration  of  1  sec. 
of  data.   The  resulting  graphs  of  this  data  were  plotted  with 
the  aid  of  the  CAL  COMP  plotter  and  a  redrawn  overlay  showing 
curves  for  selected  fan  combinations  is  given  in  Figures  25, 
25a,  26.   The  amplitudes  shown  are  relative  and  placed  on  each 
overlay  such  that  the  mean  values  of  individual  curves  are  a 
fixed  distance  (2")  apart;  these  mean  values  are  shown  on  the 
overlays  as  'tic  marks'  on  the  ordinate  axis.   It  will  be  noted 
that  the  5  fan  combination  is  missing;  during  5  fan  data  taking, 
a  slow  undulation  of  the  direct  path  sound  signal  was  noted; 
this  undulation  was  not  present  in  the  1  to  3  fan  cases  but 
was  also  noted  in  the  4  fan  case  after  about  10  minutes  at  this 
fan  combination.   It  is  surmised  that  this  variation  in  the 
direct  path  sound  was  due  either  to  relative  motion  of  the 
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source/receiver  induced  by  wave  slap  on  the  mounting  rods  or 
by  some  type  of  oscillating  current  within  the  wave  tank  itself 
For  these  reasons,  5  fan  data  were  no  longer  taken. 

The  examination  of  the  first  data  set  indicated  that 
the  time  resolution  in  the  upper  frequencies  might  not  be  fine 
enough;  consequently  an  additional  set  of  data  with  a  10  msec, 
time  separation  was  taken  encompassing  a  total  time  of  0.75 
seconds.   These  data  were  plotted  as  before  and  overlays  shown 
in  Figures  27 ,  28,  and  29. 

It  was  also  felt  that  the  total  time  display  of  data 
might  not  be  long  enough  to  show  changes  in  the  SPL  which  were 
taking  place;  therefore  a  last  set  of  experimental  data  was 
taken,  utilizing  a  20  msec,  pulse  separation  but  over  nearly 
100  points  for  a  total  time  history  on  the  order  of  2  seconds. 
These  data  were  plotted  and  overlays  shown  in  Figures  30,  31, 
and  32. 

Analysis  of  the  entire  data  array  is  contained  in  sub- 
paragraph 5  of  this  section.   Individual  plots  showing  actual 
sound  pressure  levels  and  accurate  time  scales  are  available 
in  the  Ocean  Physics  Laboratory,  Naval  Postgraduate  School. 
4 .   Post  Experimental  Data  Processing 

After  completion  of  the  graphical  presentation  of  the 
three  varied  data  sets,  and  after  being  faced  with  the  some- 
what overwhelming  mass  of  this  graphical  data,  it  was  felt  that 
further  data  processing  might  lead  to  some  simplification  of 
analysis.   The  form  that  this  further  processing  was  to  take 
was  that  of  cross-correlation;  a  method  was  desired  by  which 
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a  SPL  vs.  time  curve  at  one  frequency  could  be  compared  with 
that  of  another  frequency  and  reduced  to  a  single  number 
representing  the  relative  similarity  of  the  two  curves.  Thus, 
instead  of  correlating  in  time  or  space,  a  correlation  was  em- 
ployed using  the  algorithm  below;  it  should  be  noted  that  this 
is  a  rather  special  type  of  cross-correlation  wherein  products 
of  pressures  at  one  frequency  are  compared  with  pressures  at 

another  frequency  but  summed  over  time. 

T2 

£L        (P        (i)     -    yF       )  (P        (i)     -    yF9) 
i=Tx  Fl  1  F2  2 

Corr(F,Fn)    = 


1    2  T  T 


fl<*F  <i>  -  ^i>2  ZL 


2      h 

(PF    (i)    -   vF2)z 

i=T1         *  i=T1  2 


The  parameters  of  this  algorithm  are  defined  as  follows: 

F,   =   frequency  for  a  given  correlation  curve 

F9   =   frequency  to  which  F..  is  being  compared  or  correlated 

P   (i)  is   the  sound  pressure  (in  microbars)  at  the  ith  time 

1  interval  for  frequency  F, 

P   (i)  is  the  sound  pressure  (in  microbars)  at  the  ith  time 

2  interval  for  frequency  F2 

uF,  =  mean  value  of  P„   over  the  total  time  interval  of  concern 

1  F, 

uF„  =  mean  value  of  P_   over  tht  total  time  interval  of  concern 

2  F2 

T,   =   starting  point  in  time 
Tp   =   ending  point  in  time. 

Each  correlation  curve  for  a  given  P.  consists  of  a  set  of  11 
correlation  values  (30  kHz  to  130  kHz,  taken  at  10  kHz  inter- 
vals) .   It  should  be  noted  that  the  curves  are  normalized, 
i.e.,  when  P,  =  F2,  Corr(F1F2)  =  1. 
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The  actual  processing  of  the  data  was  relatively  simple 
since  all  data  cards  had  been  previously  punched.   A  new  com- 
puter program  was  developed  for  the  IBM- 360  to  carry  out  the 
correlation  algorithm  and  provide  the  appropriate  output.   A 
copy  of  this  program  is  contained  in  Appendix  A. 

In  addition  to  the  correlating  over  the  entire  number 
of  time  points  available,  correlation  was  also  performed  over 
the  halves  or  over  individual  fourths  of  some  curves  for  com- 
parison purposes.   Table  V  describes  those  correlations  which 
were  computed  and  plotted.   Complete  presentation  of  all  the 
correlation  curves  produced  (11  curves  for  each  case  listed 
in  Table  V)  is  not  considered  feasible  just  from  a  numerical 
standpoint  but  representative  results  are  labeled  and  shown  in 
Figures  3  3  to  39. 

5.   Analysis  of  Results 

Examination  of  Figures  25-31  gives  insight  into  the 
qualitative  aspects  of  the  effects  of  frequency  change  on  the 
scattering  of  sound.   The  relative  periodicity,  corresponding 
to  the  periodicity  of  the  wave  surface,  is  present  for  the 
relatively  low  roughness,  smoother  surface  cases.   As  the  fre- 
quency, and  hence  the  roughness,  increases  this  periodicity 
disappears.   It  is  assumed  that  the  incoherent  component  of 
the  sound  is  growing  more  predominant  and  is  responsible  for 
this  effect.   Actual  measurements  to  separate  coherent  and 
incoherent  sound  components  were  not  made.   For  the  1  and  2 
fan  overlays,  large  peaks  and  valleys  in  the  sound  amplitude 
can  be  monitored  as  they  progress  up  the  frequency  scale.   For 
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Table  V 
Correlation  Curves  Plotted 
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the  higher  fan  combinations,  the  incoherency  is  controlling, 
and  in  most  instances,  particular  peaks  and  valleys  cannot 
be  followed  over  any  appreciable  frequency  range. 

After  an  exhaustive  examination  of  the  correlation 
curves,  it  was  concluded  that  no  information  of  a  significant 
nature  could  be  discerned;  it  is  interesting  to  note  that  in 
the  3  fan  case  a  rather  rapid  de-correlation  occurs.   The  2 
fan  case,  in  contrast,  has  a  rather  high  correlation  over  the 
entire  frequency  range.   This  is  more  or  less  confirmed  by 
visually  correlating  the  curves  of  the  overlays.   It  should  be 
mentioned,  that  in  one  instance,  the  1  fan  combination  showed 
slightly  less  correlation  than  that  of  2  fans .   Also  of  inter- 
est is  the  fact  that  the  resulting  degree  of  correlation  is 
somewhat  time- independent  as  evidenced  by  the  large  degree  of 
similarity  between  the  correlations  over  the  full  number  of 
data  points  and  the  correlations  over  half  of  the  data  points. 

In  an  attempt  to  quantify  the  actual  effect  of  fre- 
quency on  the  scattered  sound,  graphs  of  the  large,  well  de- 
fined minima  taken  from  the  overlays  of  Figures  25  to  31  vs. 
time  were  plotted.   The  minima  values  were  chosen  because,  in 
general,  they  appeared  to  be  more  highly  defined  and  localized 
than  the  rather  broad  maxima  of  the  curves.   The  resultant 
minima  plots  for  the  various  overlays  appear  in  Figures  40  to 
46  and  are  appropriately  identified.   Figure  46  illustrates 
how  the  increasing  incoherency  makes  it  impossible  to  trace 
a  logical  progression  of  features  as  frequency  increases.   By 
using  an  assumed  extension  of  the  trends  developed  for  the  1 
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and  2  fan  cases,  and  also  by  assuming  that  the  maxima  are,  on 
the  average,  located  midway  between  the  plotted  minima,  it 
was  possible  to  predict  the  change  in  frequency  required  to 
go  from  a  minimum  to  a  maximum  at  any  given  instant  of  time. 
For  each  of  the  two  curves,  this   Af  was  tabulated  and  aver- 
aged.  For  the  3  and  4  fan  cases,  the  lack  of  predictability 
required  a  more  numerical  approach.   The  overlays  of  Figures 
26  and  29  (3  fans)  and  the  individual  graphs  for  the  4  fan 
cases  (not  included)  were  examined  and  the  average  frequency 
change  from  max.  to  min.  determined  and  tabulated  over  a  num- 
ber of  representative  cases.   The  results  of  both  these 
methods  are  shown  graphically  in  Figure  47.   This  curve  shows 
the  predicted  change  in  frequency  to  go  from  a  minimum  to  a 
maximum  plotted  against  the  average  wave  height  for  a  given 
fan  combination.   The  1  and  2  fan  frequency  changes  were 

calculated  from  10  values  and  the  3  and  4  fan  cases  from  14 
and  17  values  respectively.   Standard  deviations  about  each 
of  these  points  are  also  shown  in  Figure  47.   The  curve  appears 
linear  for  the  smoother  surfaces  progressing  into  non-linearity 
by  the  time  the  4  fan  surface  heights  are  reached;  however 
this  lack  of  linearity  of  the  last  point  may  simply  be  due  to 
a  lack  of  resolution  below  10  kHz  increments . 
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V.   FUTURE  RESEARCH 

The  author  feels  that  the  research  described  in  the  fore- 
going pages  is  a  viable  means  for  gaining  new  knowledge  and 
insight  into  the  basic  process  of  scattering  of  sound  off  of 
a  statistically  rough  surface.   Further  investigation  should 
provide  additional  information;  a  number  of  approaches  to 
improve  the  previous  research  are  recommended  below: 

1.  Investigation  into  the  feasibility  of  using  a  5  kHz 
base  frequency,  instead  of  the  10  kHz  used  herein,  should  be 
conducted.   A  few  trials  along  these  lines  were  attempted 
with  indications  that  the  method  is  feasible.   The  frequency 
resolution  would  then  be  twice  that  given  in  this  paper  and 
should  be  particularly  valuable  in  giving  a  closer  look  at 
the  3  and  4  fan  surfaces  at  lower  roughness  values.   Rough- 
ness values  less  than  1.0  should  specifically  be  looked  at. 

2.  A  method  for  separately  looking  at  the  coherent  and 
incoherent  components  of  the  scattered  sound  should  be  de- 
veloped.  This  appears  to  be  capable  of  being  accomplished 
using  the  phase  information  present  in  the  FFT  algorithm  used 
in  the  digital  computer.   The  SKIP-1C  program  has  this  output 
available.   A  possibility  of  merging  the  capabilities  and 
outputs  of  the  SKIP-1C  and  THCDB(l)  programs  should  be 
attempted.   By  removing  the  incoherent  component  at  the  higher 
surface  roughness  values,  more  information  can  be  attained 
from  the  chaos  which  currently  presents  itself. 
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3.  Alternate  means  of  presenting  correlation  data  should 
be  investigated;  one  possibility  would  be  to  calculate  the 
correlation  of  excess  pressure  amplitudes  relative  to  the  mean 
value  for  a  given  frequency  increment  at  a  given  instant  in 
time  and  to  display  the  resultant  correlation  values  as  a 
function  of  time.   The  advantage  of  this  method  is  that  it 
would  lend  further  support  to  the  hypothesis  of  consistent 
frequency  separation  (delta  f)  between  relative  maxima  and 
minima  discussed  in  section  IV. B. 5.   It  is  anticipated  that 
such  curves  would  reveal  a  large  negative  correlation  when  the 
appropriate  frequency  of  Figure  47  is  reached.   Correlation 

at  twice  this  frequency  increment  should  result  in  large  posi- 
tive correlation  values  once  again.   Analysis  of  existing 
data  arrays  in  this  manner  is  feasible  due  to  the  prepunched 
data  cards  available  in  the  Ocean  Physics  Laboratory. 

4.  The  effects  of  variation  of  geometry  to  include 
larger  angles  of  incidence  should  be  investigated. 
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APPENDIX  A:   PROGRAM  DESCRIPTIONS 


SKIP-1C: 

Input: 

1 

2 

3 

4 

Options : 

1 

2 

Output: 

1 

3. 


Sample  Frequency 
Number  of  Sample  Points  per  FFT 
Number  of  FFT's 

Band  pass  limits  desired  on  printout 
DC  filtering 

Triangular  hanning  function 
Individual  power  spectral  densities 
(PSD)  for  each  FFT 

Normalized  average  PSD  vs.  frequency 
over  the  stipulated  bandpass 
Average  phase  vs.  frequency  over  the 
stipulated  bandpass. 
For  each  FFT,  the  maximum  PSD  and 
its  associated  frequency. 


CAL-01: 


Input: 


Output 


1.  Number  of  samples  (maximum  819  2) 

2.  Input  reference  voltage  set  at 
1.000  volts  (RMS) 

3.  Voltage  signal  from  wavetank  probe 
1.   VOUT  =  AC  voltage  of  wavetank  probe 
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THCDB: 


Input: 


1.  Number  of  sample  points  per  FFT 

2.  Total  number  of  FFT's  desired  (maximum 
of  50) 

3.  Conversion  factor  for  transforming 
PSD  to  volts 

4.  Amplifier  gain 

5.  Line  loss 

6.  Hydrophone  sensitivities  for  every 
5  kHz  from  0  -  160  kHz 


Options 


Output : 


THCDB (1)  : 


1.   SPL  or  VL  output 
Noise  run 

1.  Mean  sound  pressure  level  (or  voltage 
level)  vs.  frequency  (5  kHz  increments 
from  0  -  160  kHz) 

2.  Standard  deviation  of  SPL  vs.  frequency 

3.  Signal  to  noise  ratio  vs.  frequency,  if 
previous  noise  run  has  been  performed 

4 .  Output  can  be  stored  on  a  digital 
cassette  or  printed  on  the  terminal 
printer 

Same  as  the  previous  program  except  the 
program  will  accept  raw  A/D  input  from  a 
digital  cassette  and  it  can  handle  up  to 
100  FFT's  as  the  associated  output. 
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PERK-1A: 


Input 


Output : 


PRK-MOD: 


Program  to  perform  A/D  sampling,  delay  a 
specified  time  and  write  the  A/D  data  on 
a  digital  cassette  (to  be  used  as  input 
for  THCDB(l) ) . 

1.  Total  number  of  A/D  samples  (maximum  of 
8200) 

2.  Number  of  samples  per  FFT 

3.  Sample  frequency 

4.  Number  of  FFT's  desired 

5.  Desired  time  delay  between  FFT's 

1.   Raw  A/D  data  stored  on  digital  cassette. 

Same  as  previous  program  except  total 
number  of  points  has  been  increased  to 
a  maximum  value  of  13,000 


SO  RTF : 


Program  to  take  the  output  of  THCDB(l) 
and  sort  the  resulting  FFT  data  by  fre- 
quency.  It  can  accept  up  to  100  blocks 
of  FFT  data  with  up  to  75  frequency 
values  in  each  block. 

Output  is  SPL  vs.  time  for  each  frequency 
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